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Abstract. In numerical experiments involving nonlinear solitary waves propagat- 
ing through nonhomogcneous media one observes "breathing" in the sense of the 
amplitude of the wave going up and down on a much faster scale than the motion 
of the wave - see Fig. [5]below. In this paper we investigate this phenomenon in the 
simplest case of stationary waves in which the evolution corresponds to relaxation 
to a nonlinear ground state. The particular model is the popular So impurity in 
the cubic nonlinear Schrodinger equation on the line. We give asymptotics of the 
amplitude on a finite but relevant time interval and show their remarkable agree- 
ment with numerical experiments, see Fig. [1] We stress the nonlinear origin of the 
"breathing patterns" caused by the selection of the ground state depending on the 
initial data, and by the non-normality of the linearized operator. 



We study a simple model of relaxation to a nonlinear ground state. Our equation is 
the one dimensional nonlinear cubic Schrodinger equation with a small delta potential: 



where < |g| <^ 1. The nonlinear ground state minimizes the corresponding energy 
(12. ip for a prescribed norm, and is explicitly given by 

(1.2) vx{x) = Xsech{X\x\ +tcinh'^{q/X)) , \\vx\\l = 2{X - q) , A > |g| . 

A simple rescaling allows the reduction to the case A = 1 and we obtain 

Theorem 1. Suppose thatu{x,t) solves (II. ip . u{x,0) G H^{M.) is real and even, and 
that 



1. Introduction 







where w{x,t) is given explicitly in (15. 6p . In particular, for 1 <^t < C\q\ '^l'' 
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Figure 1. Breathing patterns for u{x,0) = sech(a;/(l + g))/(l + q) 
(the initial data is rescaled so that the ground state to which it relaxes 
is Vi{x)): the plots show |M(0,t)| and the asymptotic prediction (11. 5p 
given in Theorem [H for q = 0.05 and q = 0.01. The agreement is 
remarkably good for times much longer than given in the theoretical 
result. 



The conditions on wq in (II. 3p can be weakened considerably, and in particular we 
only need estimates for k ,£ < N for some A^. Theorem [2] below gives a statement 
which depends only on ||wo||hi being small and the more explicit results in Theorem 
[Hcome from our close analysis of the propagator exp{—itCq^\) appearing in (13. 2p . As 
explained below our motivation comes from the study of solitons and the intial data 
in which we are most interested is u{x, 0) = sechx. For that case, the comparison of 
the theorem with numerical results is shown in FiglH 

The reduction to the case A = 1, mentioned before the statement of Theorem [1] is 
straightforward: let u{x,t) = A~^m(A"^x, A^^t). Then u solves (ll.ip with q replaced 
by \"^q, idtu + ^d^u + qX~^6o{x)u + ulul"^ = 0. Now, if we suppose the theorem holds 
when applied to u replacing u and q/X replacing q, then we can deduce the theorem 
in its current form. Thus, it suffices to prove the A = 1 case. 
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In the remainder of the introduction we will discuss our motivation and relations to 
existing literature, possible approaches to obtaining finer asymptotics, and a simple 
example of a breathing pattern for nonnormal operators. 

1.1. Motivation. Mathematical studies of relaxation to ground states for nonlinear 
Schrodinger equations have been recently conducted in a number of mathematical pa- 
pers, see Soffer-Weinstein [26], Tsai-Yau [28], Gang-Sigal [H], Gang-Weinstein [T2] . 
and references given there. The particular focus is on the behaviour as t ^ oo (gen- 
uine relaxation in the sense of pure mathematics) and the allowed non-linearities typ- 
ically exclude standard examples from the physical literature. For the cubic nonlinear 
Schrodinger equation (NLS) on the line, that is for (11. ip with q = 0, the nonlinear 
relaxation can be studied in great detail using methods of inverse scattering theory 
pioneered by Zakharov-Shabat [31] - see Deift-Its-Zhou [7], Deift-Zhou [6], and [HJ 
Appendix B] for recent advances and references. The case of g 7^ with even initial 
data is also in principle accessible by these methods as was pointed out by Fokas [9]. 

The goals of this paper are more modest: we explain a phenomeno logical fact 
occuring on shorter time scales for a simple physically relevant model of NLS with 
small 6 impurities. Numerous references for this model in the physics literature can 
for instance be found in [1] (where it is used to model more realistic narrow traps), 
[5],[l3], and [22]. See also [21] for a recent numerical study and further pointers to 
the literature. 

Our motivation came from observing a common phenomenon illustrated in FigEl 
In [16] we have shown that the solution of ( 11. ip with u{x, 0) = sech(x — ao)e^^'"° 
satisfies 

(1.6) \Ht, .) - e^(— W)-We^7{t) _ amimm < C\q\'-'' , 

for < t < S{vq + |g|)~^/^ log(l/|g|), < |g| <^ 1, and where a, v, and 7 solve the 
following system of equations 

^a = v, ^v = ]-qd^ (sech^) (a) , 
at at I 

^ ' d I v'^ 1 

— 7 = - + y + gsech^(a) + -qad^{sech.^){a) , 

with initial data (oq, Vq, 0) (please note that the sign convention for q has been changed 
here). As was pointed out there, as seen from explicit constants in coercive estimates, 
these asymptotics require q < 0.01 to hold accurately. From the semiclassical point 
of view q = K^, where h is the effective Planck constant of the problem, so that means 
/i < 0.1 - see [T7] for an explanation of this scaling philosophy. 

In Fig. [2] the dashed line shows the motion of the center of the soliton in the case of 
q = 0.05 (which is a borderline case for the applicability of (11. 6p ). We see oscillations 
with the period proportional to g"^/^ in agreement with (11.70 . The continuous line 
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Figure 2. Fast and slow oscillations in the motion of the soliton with 
initial condition sech(x+3) moving in the field of —5q{x)/2Q (top figure) 
and — (5o(a^)/40 (bottom figure). The center of the soliton oscillates 
around a; = with a much larger period than the rescaled amplitude 
A{t) = 30{\u{a{t),t) \ — 1), where a{t) is the center of the moving of the 
moving soliton. The periods of the fast oscillations are close. 

shows the oscillation of the amplitude: we look at the deviations of the value of the 
solution at the maximum of |u(x,t)| in x from 1, the maximal value of the absolute 
value of the soliton solution. The oscillations are much faster than the oscillations of 
the center of the soliton and the period is close to being fixed. Numerical observations 
suggest that the period is almost independent of q. 

This "breathing" behaviour is even more striking in movies of numerical solutions 



see for instance the last movie in http: //math.berkeley . edu/~zworski/msg.pdf ). 



The slowing down of the soliton and the shedding of its mass seem closely related to 
these breathing patterns. 

As the first step to understand solitons moving in nonhomogeneous media we study 
the stationary case, that is (11.11) with initial data given by u{x,0) = sechx. The 
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results of [16] recalled in (11. 6p and (11.71) show that 

.) _ enW sech(.) ||hi(e) < , 

for < t < log(l/|g|), and ■y(t) = (1/2 + q)t. In fact, an application of the 

method of [16] shows that for some 7(t), 

||n(t,.)-e^^Wsech(.)||Hi(M) <C|g|, 

for all times. Here we could replace sech with vx for any A = 1 + 0{q). 

Hence the breathing patterns must involve higher order asymptotics and since 
|g| = h'^, the natural next step is \q\'^^'^ = h?. Theorem [1] provides that next step 
on a time scale which allows seeing a large number of oscillations. The numerical 
experiments show a very good agreement with asymptotics provided by (II. 5p and 
suggest that they are valid for times longer than t <C 

Finer asymptotics might be possible if one adapts some of the methods of [26], |llj . 
and [12], but it is not clear which direction should be taken for the efficient study of 
moving solitons. We opted for the simplest at this early stage. 

1.2. Nonlinear aspects of "breathing". We first compare the breathing patterns 
observed here with amplitude oscillations in the relaxation to the ground state of a 
linear problem, iut = —Uxxf^ — qSqU, < g -C 1. If the initial data is equal to Uo{x) 
and is real and even, a heuristic approximation for the solution is 



Vi • 

Although the zero resonance disappears (see [23] and §4.ip . for q small we expect the 
behaviour 1 / to persist for long times - see §5.4[ 

This is very different from (11.51) . The main difference is that in the nonlinear prob- 
lem the eigenvalue is not fixed but it is selected depending on the initial condition. 
The approximate selection is given by the formula for A in Theorem [1] and a more 
precise selection method is given in Proposition 13 . II preceding Theorem O In particu- 
lar, the periods oscillation for a fixed initial condition are approximately independent 
of q. Since the linear eigenvalue is fixed and depends on q this is strikingly different 
in linear relaxation - see FiglH 

The origin of the phase in the second term in (II. 5p lies in the properties of the 
non- normal linearized operator for (II. ip . and in particular in the coupling responsible 
for the nonnormality. We do not yet have a fully conceptual explanation for that 
other than the analysis of (15.60 . 

We present a simple example illustrating how non-normality can be responsible 
for "breathing", that is oscillations in the amplitute, absent for normal operators. 
Suppose that G M (note that we allow both positive and negative a, f3), w = 
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Figure 3. Examples of linear relaxation: the initial condition is 
u{x,0) = secha; and the potentials are —q6Q{x) withg = 1/2 for the 
top graph and q = for the bottom graph. We expect the pe- 

riods of oscillations to be Air/q'^. Consequently, changing g = 1/2 to 
q/\/2 = \/2/2 we expect the period to double and to see that we plot 
\u{0, t) \ in the top graph and |m(0, t/2) \ in the bottom graph: the agree- 
ment of the periods is striking. The horizontal lines correspond to the 
asymptotic values Vq{0) J sech.xvq{x)dx, Vq = y/qexp{—q\x\) . 



[Rew Imw] 



'1.8) R 



and we consider the evolution 



a + dl 



;i.9) 



dfW = Rw 



We consider R as an operator on iir^(M; C) x C) and write out explicit formulas 

for the complex-valued vector plane-wave solutions associated to (generalized) eigen- 
values iicj of R. From these, we build real- valued plane- wave solutions w = [wi ^2]^ 
to the matrix equation (11.91) . When these are converted to complex numbers as 
w = Wi + iw2, we find that if a = (3, then w is unimodular but if a 7^ /3, then w is 
not unimodular and we see oscillations in amplitude. 

Let uj > 0. We seek complex- valued (vector) plane- wave solutions to (11.91) with 
generalized eigenvalue ±iuj. Let 7 > ma.x{^ymaJ<.{a, 0), v/max(/3, 0)) > be the 
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unique solution to 



Set 



Now let 



1 

ia 



a 



a 



72-/3 



5(7) 



1 

-zcr 



Then f (7), f (— 7) are two plane- wave solutions to Rv = iuov and f (7), f (— 7) are two 
plane-wave solutions to Rv = —iuov. From this we see that 



t<±7), 



w(±7) 



are four solutions to (11. 9p . Thus 

cos(tCi; + 7x) 



w 



-7) 



2^ 



a sin(tci; + 7a;) 

is a real solution. Forming a complex number from this vector, we obtain 

cos(tco' + 7x) + ia sin(tu; + 7x) 

which has constant-in-time modulus if and only if cr = 1. 

The linearization of (11.11) around the solution vi gives the non-normal operator 
defined below in (I2.12p . The motivation for studying the operator R above is that 
has the form of R with a = (3 = —1 for |x| large but with a = 5, /? = 1 for |x| near 0. 

1.3. Method of proof and organization of the paper. The proof Theorem [1] 
consists of two parts. The first is a nonlinear perturbation theory presented in Theo- 
rem [2] and gives an approximation of the solution by the linearized flow. The second 
part is the precise analysis of that linearized flow on time scales consistent with the 
approximation given in Theorem [21 

In ^we present various standard facts about the nonlinear Schrodinger flow with 
an external delta function potential. The Hamiltonian structure of this flow with 
respect to the symplectic form uj{u,v) = Im J uv is particularly crucial. It plays an 
important role in ^ where it is used to select the nonlinear eigenvalue of the limiting 
"relaxed state" . The other component of the proof of Theorem |2l which is the main 
result of that section, is the coercivity estimate allowing the control of norm by 
the linearized operator. The estimates on the propagator are similar to the estimates 
in §5 of [T7j: instead of the L^-energy method we estimate dt{CqU,u), where Cg is 
essentially the Hessian of the Hamiltonian (see (12.81) ). The initial data is assumed to 
be even which is particularly important in the case of g < (repulsive 6 potential) as 
the ground state is then unstable - see [22] . 
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In ^we recall Kaup's explicit spectral decomposition of the linearized operator for 
the focusing cubic NLS on the line. As shown in Appendix [A] that basis can also be 
discovered via simple numerical experimentation. We use it to obtain a representation 
of the propagator and apply it to see the relaxation to solitons in the free case. For 
initial data close to the solitons this crude approximation is remarkably close to the 
precise results given by full inverse spectral method - see FigJSl 

The results of §l]lead to an almost explicit spectral decomposition for the operator 
with < |g| ^ 1 - that is presented in ^ and Appendices [B] and O We follow 
the general theory of Buslaev-Perelman and Krieger-Schlag in that particular setting. 
More general non-linearities could also be allowed but since we are ultimately inter- 
ested in the comparison with numerics the explicit nature of Kaup's basis is very 
useful. As in §1] this leads to a representation of the propagator, (15.61) . used in the 
statement of Theorem [H An asymptotic analysis of §5.41 gives the approximation 
f fT3D illustrated in FiglH 

Acknowledgments. We would like to thank Michael Weinstein and Gahna Perel- 
man for stimulating conversations and e-mail exchanges. The work of the first author 
was supported in part by an NSF postdoctoral fellowship, and that of the second 
second author by the NSF grant DMS-0654436. 

2. Preliminaries 

In this section we review various basic aspects of the equation (II. ip . 

2.1. Hamiltonian structure. The nonlinear Schrodinger equation (11.11) describes 
the Hamiltonian flow on on H^(R, C) for the Hamiltonian 

(2.1) H,{v)'^\J{\d.v\'-\v\')dx-\q\vm\ 
More precisely, we consider 

as a real Hilbert space with the inner product and the symplectic form given by 

(2.2) (m, =^Re j uv , uj{u,v) {iu,v) = Ira 

Let Hq given by (12.11) . or be a more general function, H : V ^ R. The associated 
Hamiltonian vector field is a map '■ V — > TV, which means that for a particular 
point u E V, we have {'Eh)u G TuV. The vector field Eh is defined by the relation 

(2.3) iuiv, iEH)u) = duH{v) , 
where v G TuV, and duH : T^V ^ M is defined by 



duH{v) = ^ 
ds 



Hiu + sv) . 

=0 
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In the notation above 

(2.4) dH^iv) = {dH^,v) , {EH)u = -dH^. 



For H = Hq given by fl2.1l) we compute 

duH{v) = Re y {{l/2)d^ud.A' - \u\'^uv)dx - Re{qu{0)v) 

Thus, in view of ([23D and Q, 

The flow associated to this vector field (Hamiltonian flow) is 



-d^u — \u\ M — q6Q{x)u 



(2.5) it = (Eh)u = ^ y^^^l^ ~ 1^1^" ~ qSo{x)uj . 

2.2. Well posedness in H^. The discussion here has been formal but it is well known 
that the equation (11.11) has global solutions in for more general nonlinearities, 
Iwp^^M, 1 < j9 < 5. For the reader's convenience we recall the standard argument. 
We have the following basic estimates: 

(2.6) 11^11^00 < C||M||i2||-u'||i2 , 

(which follows from the fundamental theorem of calculus: m(x)^ = 2u{y)u'{y)dy) 
and thus 

1 p-l p+3 1 2(p+3) 

J- n_1 II iiO , ^11 / II ^ II II ^ , -L II /II O ^/ 11 11 ~E ~ 



\U\ 



^+'< ^\\U\\1-J\\U\\1,<C\\U'\\J \\U\\A < -:\\u'\\^r2+C'\\uLr'' 



|KO)r<^||n'||i. + Cg^|Hli.. 



Hence, 



H,iu) = \\\urL'^-^^\\uci-^\um' 

1 2(p+3) 
-i-ii /iio ^11 II — E — 



and consequently. 



>-\\u'\\l-C\\u\\^l^ -Cq'\\u\\l, 



2(p+3) 

\u\\li < 8Hg{u) + CWuW^l-" + (Cq^ + l)||M||i2 . 



Since the energy, Hq{u), and mass, ||«||^2, are conserved, we see that if the solution 
exists in H^, its norm is uniformly bounded. Thus we only need to show local 



10 J. HOLMER AND M. ZWORSKI 

existence in H^. Let us fix T > and, for u = u{x,t), define the norm 

sup \\u{»,t)\\H^ . 

0<t<T 

Solving fll.ll) . u{x, 0) = uq{x), is equivalent to finding the fixed point of the operator 

« Jo 

Here the operator exp{it{dl/2 + q6Q{x))) is unitary on (see the discussion of the 
operator L given in (I2.9P below) and preserves 

H,{n) = \\\u'\\l^-l\um\ 

Again using (12. 6p . 

IWu'Wl - Cq'Ml, < HMt)) < + Cq'Ml, . 

Therefore, if u{t) = exp(it(9^/2 + g5o(x)))Mo, 

lh'ml.<H,{u{t)) + Cq'Mt)\\l, 
= Hg{uo) + Cg^lluollia 
< ^ll^olli^ + Cq^\\uo\\l2 

From this, we see that 

is bounded with norm independent of t. This and the estimate 

\\\u\P~^U-\vf~'\\Hi < C{\\\uf^'\\H^ + \\\v\''~'\\m)\\u-v\\Hl 

< C{\\u\\hi + \\v\\m)P~^\\u-v\\H^ . 

give 

||$(u) - $(i;)|U < CT{\\u\\x + \\v\\xY ^\\\u - v\\x) 
so that for T small fixed point arguments can be used to obtain a solution in H^. 

2.3. Nonlinear ground states. The minimizers with a prescribed norm are given 
by critical points of the Hamiltonian with the constraint added (and playing the 
role of the Lagrange multiplier): 

(2.7) SM = H,{u) + ^\\u\\l,. 

Then for the ground state given by (II. 2p we obtain 
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where the Hessian is the following self-adjoint operator on iJ^(M, C) ~ iJ^(R, R) © 

" 

V. ' 

Qv^ - 2q5o) , 
2v^ - 2q5o) . 

In view of the 6 functions, the definition of the operators Lg± is given by choosing 
the correct domain for the operator. To see what it is let us first examine the basic 
case of 

(2.9) L = -dl + V -q6o, 

on M, where is a smooth real- valued potential, rapidly decaying at oo. Suppose 
that u G L^(M), Lu = /, and / G L^(M). This implies that away from we have that 
dlu G and thus u G i/2(M\{0}). In order that / remain a function across x = 0, 
we must have that u(x) is continuous at x = and 

(2.10) n'(0+) - n'(O-) = -qu{0) 
Thus a natural domain to consider for L is 

(2.11) V = {u\uE H'^(R\{0}) , u is continuous at x = and (KTO\i holds} 

By verifying that the operators L ±i are both symmetric and surjective on V we see 
that L is self-adjoint with domain V. 

2.4. Linearization and the Hamiltonian map. For S : V ^ satisfying S'{u) = 
we can invariantly define the Hamiltonian map, 

using the well defined Hessian of £^ at m: 

{E"{u)X,Y) =uj{J,TX). 

In other words, the Hamiltonian map is the linearization of the Hamilton vector field 
of E. See for instance [HI Sect. 21. 5] for a general discussion, and [171 Lemmas 2.1, 
2.2] for relevant facts in our context. 

For V = (M, C) with the symplectic form (12. 2p we have 

= -i£" , 

and for £ given by ( 12. 7p we have 

(2.12) = -iC, = 



(2.8) Cq = 



where 

V = ^(A^-5^ 

V = ^(A^-5^ 





-V 
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The matrix representation is based on the identification 

It is also convenient to consider the equivalent matrix representation using the iden- 
tification, 

if ^(M, C)3u - [u, m]* G A C H\R, Cf , 

which gives 



u 



dcf 



1 



V2 



(2.13) 



H, 



q,X 



-dl + 



dl + A2 







dl- 



dl - \^ 





'-4 -2" 


+ sech^ X 




2 4_ 



-4 
2 



(when there is no, or little, chance of confusion we supress A in our notation; most of 
the time its value is taken to be 1). This representation is convenient when we study 
the spectral decomposition of Fg. The factor | was introduced to make the notation 
simpler and to have a better agreement with the standard notation of [TH], P, [2T]. 

Since the energy Hq{u) differs from 2Sq^x by the additive mass term, A^||u|p/2, 
these linearizations differ from the linearization of fll.ll) by a constant only. 

For future reference we also note the symmetries of Hq. Let aj be the Pauli matrices. 



(2.14) 



dcf 


'0 


1" 


def 


'0 -i 


dcf 


'1 " 


ai = 


1 





, 0-2 = 


i 




-1 



We recall that they are characterized by the properties that (t| 
Using this notation, 

(2.15) cTiFgdi = -H, 



I and a* = aj. 



a-3if„a3 = H* 



General considerations show that cr{Hq) C M U iM, and in fact cr{Hq) \ Cpp (Hq) = 
(— oo, —1] U [1, +oo). The fact that all pure point spectrum is contained in M U 
follows by examining the squared operator which turns out to be self-adjoint, (see 
Buslaev-Perelman [H §2.2.3]). This can be done despite the fact that the operator 
contains the 6o potential. To see that consider again the operator L given in (12.91) 
and suppose that we want to consider L^, the squared operator. Away from x = 
0, we see that we must have u G if^(M\{0}). If Lu = /, then we need f E T). 
Since / is continuous at 0, we see from the equation Lu = f that lim^^o- d^u^x) = 
\im:c^o+ dlu{x) . Moreover, taking 



u"iO) 



dcf 



lim dlu{x) 



lim dlu(x) 
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implies 

u'\Q) = V{0)uiO) - /(O) 
Away from x = 0, we have —u'" + V'u + Vu' = f and thus the condition that 
/'(0+) - /'(0-) = -g/(0) becomes 

(2.16) u"{0+) - u"{0-) = -qu"{0) 

Define 
(2.17) 

V = {u\ue i/^(R\{0}) , u, u" are continuous at x = and fl230|) . fl236|) hold} 

Provided u E V, LF'u is defined and belongs to L^(]R) - so there is no need to worry 
about the square of the delta function not being defined. Indeed, as soon as we know 
that M e P as defined here, then one need only compute {—d1 + VYu away from x = 
to obtain LP'u. It is thus natural to consider the squared operator onVxV, and 
Hg can in fact be shown to be self-adjoint on this domain. 

2.5. Symmetries and the generalized kernel. As in [16] and [17] it is convenient 
to introduce a natural group action on H^: 

^2 H'3u^g-ueH\ (g ■ u){x) e'^ e''^'^'-^^ fiu{f,{x - a)) , 

g = {a,v,'y,fi) G x IR+ . 

This action gives a group structure on M.^ x ]R_|_ and it is easy to check that this 
transformation group is a semidirect product of the Heisenberg group and M+: 

G = H3>^R+, /i ■ (a,w,7) = (-,/i?;,7) . 

The Lie algebra of G, denoted by g, is generated by 61,62,63,64, which in the 
infinitesimal representation obtained from f l2.18p is given by 

(2.19) 61 = -^2;, 62 = ix, 63 = 2, ei = d^-x. 

It acts, for instance, on iS(M) C H^, and by X G g we will denote a linear combination 
of the operators Cj. We note that for g = (and hence v = vx = X sech(Ax)) 

uj(ei ■ 62 ■ w) = 1 , uj(e3 ■ 64 ■ f ) = 1 , 

(2.20) 

uj{ej ■ 63 ■ f ) = uj{ej ■ 64 ■ t") = , j = 1, 2 . 

In the case of g = the Hamilton vector fields, H/^-g, ^ = + A/4, are tangent to 
the manifold of solitons, G ■ v - see [lOl §3] or [TTJ §2.2]. Hence JFq = —iCo preserves 
Ty{G ■ w) ^ g ■ t>. In fact, g ■ f is the generalized kernel of —iCo: 

iCo(ei ■ v) = , ijCo(e2 ■ v) = ei ■ v , 

(2.21) 

iCo{e3-v) = 0, iCo{e^ ■ v) = 63 ■ V . 

The first and third equation are an immediate consequence of the invariance of solu- 
tions under the circle action u 1— > e^^u, and translation in x. 
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Figure 4. The spectrum of the operator Hq. In the notation of §2.51 
the generahzed eigenspace at is spanned by Cj ■ sech, j = 1, ■ ■ ■ , 4. 

When g 7^ we lose the translation invariance but we still have iCq{ei ■ v) = due 
to the preserved circle action symmetry. As shown in Appendix [B] the generalized 
kernel is given by 



Vo = span{w3,t;4} 



where 



V3[X) = tVx[X) 



, ^^4(2;) = dx vx{x), iCgVs = , iCgV4 = V3 . 

1 X=l 

Hence vj, j = 3, 4, are the generalizations of ej ■ v, and in fact, 

f 3 = 63 ■ f , f 4 = 64 ■ w + Ci/i (q) , 

2.6. Coercivity estimate. Finally we recall the crucial coercivity estimate which 
in a more general form is well known since the work of Weinstein [29] . For the special 
case at hand an elementary presentation can be found in [161 §4]- 

For g = we have the following estimate: Let w G if^(M, C) and suppose that for 
any X E Q, uj{w, X ■ t]) = 0. Then, 

(2.22) {CqW, w)> - 0.0555||w||^i , po - ^ 



7 + 2po 2(12 + 7r2 
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3. Nonlinear perturbation theory 

In this section we prove a result describing eigenstate selection and nonlinear flow 
approximation for a time depending on the initial data and on the size of q. Although 
we restrict our attention to the physical (and completely integrable) case of the cu- 
bic NLS the arguments apply to nonlinearities for which the Weinstein coercivity 
conditions are satisfied (see [2^ and Lemma [3.41 below) . 

Recall 

vx,q{x) = Asech(A|x| + tanh"^(g/A)) , ||fA,g(a;)||i2 = 2(A - q) . 
Define the projection 

(3.1) P\,qip = Uj{ip, dxVx,q)iVx,g - Uj{(p, iVx,g)dxVx,g 

onto the generalized kernel 

Vx,q = spanB;{n'A,g, dxVx,q} 

of Cx,q- We will only use the q subscript when it is needed for clarity (e.g. in the 
scaling argument below). We will also drop the A-subscript when A = 1. Recall that 
{u, v) = Re J uv. 

Proposition 3.1 (Symplectic orthogonality). There exists 6 > such that the fol- 
lowing holds. IfipE and there exists Aq > 0, 6*0 G M such that Hy? — e'^^'fAoll//! ^ ^, 
then there exists A G (0, +oo), 6' G M such that Px{e~'^^ip — Vx) = 0. 



Proof. Let F : x (0, +oo) x R ^ R"^ he given by 

F{uo,\,9) 

Fix 6*0, Aq. Note that F(e*^"t%, Aq, 6^0) = and the matrix 

[dxF deF] 



uj[uo — e Vx, le Vx) 
uj{uo - e'^vx, e'^dxvx) 



1 
-1 



is (uniformly in A, &) nondegenerate at uq = e*^°WAo, A = Aq, 6* = Qq. The implicit 
function theorem completes the proof. □ 

Theorem 2 (Nonlinear perturbation theory). Let I d (0, +00) and |g| <^ 1. Suppose 
that u{x,t) is an even solution to (11. ip and Wo{x) '= u{x,0) — e'^^vx{x) satisfies 
\\i^o\\m <h<^l and Px{e~'^Wo) = for some Xe I, 9 eR. Then 

(3.2) \\u{t) - e**^'/2 (yx + e-**^^-wo) ll^i < Ct(l + t)h^ , 

for allO<t<^ 
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Remark 3.2. The constant C depends on / (the range of values in which A hes), and 
the restrictions |g| <^ 1, /i ^ 1 and t ^ h^^^'^ all indicate an implicit (small) constant 
depending on /. 

Remark 3.3. We will ultimately take h = Cq to prove Theorem 1 in §5.51 Note 
that our use of h here is different from the connection to the semiclassical problem 
(discussed in the introduction) where q = h?. 

Lemma 3.4 (Coercivity). There exists Cq > (independent of q) with the following 
property: If \q\ ^ 1, Pqf = and f is even, then 

11/11^1 <co(£ J,/) 

Proof. We have, with / = /i + /2 

(£J,/) = (L,+ /l,/l) + (V/2,/2), 

where Lg+ and Lq_ are the self-adjoint operators defined in (12. 8p . It suffices to prove 
that if / is even and real-valued, then 

(3.3) {f,v) = =^ (L+/,/)>c||/||i., 
and 

(3.4) {f,d,v,\,=,) = =^ {L^fJ)>c\\f\\l.. 

The operators L± (defined as Lg± with q = 0) were analysed in [16], §4], and it was 
proved there that cr(L_|_) = { — |, 0} U [|, +oo) and C3"(L_) = {0} U [|, -|-oo). Moreover, 
the eigenvalues and normalized eigenfunctions are given explicitly: 

L+(^sech2) = -K^sech^), L+{^sech.') = 0, L_(^sech) = 0. 

By perturbation theory (this is standard perturbation theory for 2nd-order scalar self- 
adjoint operators, as opposed to the perturbation theory of Appendix [B]), a{Lq^) = 
{Ai, Ao} U [|, +oo), where Ai = — | + 0{q) and Aq = 0{q). Moreover, the normal- 
ized associated eigenfunctions, gi and go, 

Lq+gi = Aif?!, Lq^go = Aof^O ! 

satisfy 

gi{x) = ^ sech^ x + 0{q), goix) = sech' x + 0{q) . 

In particular, (7o(x) is "nearly" odd. However, we also have that Lq+go{—x) = 
\ogo{—x), and hence go{—x) = cg{x) for some constant c. If go{0) ^ 0, then c = 1 
and hence g^ix) is even which contradicts the fact that it is nearly odd. From this we 
conclude (7o(0) = 0, and since g^ix) is not identically zero and solves a second order 
ODE, we must have 5^0(0) 7^ 0- Taking the derivative of the identity g^i—x^ = cgo{x) 
and evaluating at x = gives that c = — 1. Hence go is exactly odd. 
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Now we prove (13.31) . Since / is assumed even, we have by the Spectral Theorem 
that if = 0, then {Lg+fJ) > ^\\f\\l2. By [HI Lemma 4.2] (with, in the 

notation of that Lemma, vq = gi, vi = v, cq = — Ai), we have that if (/, v) = 0, then 



2 

L2 



By the perturbation theory, the coefficient evaluates to ^ — | + (9(g) > 0. 

Now we carry out the analysis of Lq_. By perturbation theory, we know that 
a{Lq^) = {A2} U [|,+cxd), where A2 = 0{q). However, direct calculation shows that 
Lq-V = (where v{x) = sech(|x| + tanh~^g)), and thus A2 = 0. By the Spectral 
Theorem, if {f,v) = 0, then (Lq^fJ) > |||/||i2. By [HI Lemma 4.2], we have that 
if {f,dxvx\x=i) = 0, 



Elliptic regularity completes the argument. □ 
Proof of Theorem\^ Let u{x,t) = e~^^X~^u{X~^x, X~'^t). Then u solves 

idtu + \dlu + ^6ou + = 
A 

( (ll.ip with q replaced by q/X) with initial data Uo{x) = e~*^A~^no(A~^x). Moreover, 
11% - e'\!x,q\\H^ <h^ \\uo - vi^q/xWm < C{X)h 

and 

Px,e,q{uo - e'%x) = ^ Pi,o,q/x{uo - Vi,g/x) = 
Hence, it suffices to prove the theorem in the case A = 1, 6' = 0. Let 

vJx) = dx 



V3[X) = lVx[X) 



Vx[X 

X=l 



X=l 

Note (by direct computation) that iCqV^ = and z£gf4 = and 

Pw = Uj{w, f4)t'3 — uj{w, ^3)^4 

is the symplectic orthogonal projection onto the generalized kernel Vq = spanjt'a, V4}. 
Define v{x) =^ '^x{^)\x=i ^^'^ '"^(^) relation u(t) = e**/^(t; + w(t)), and then 



dcf 

note that w solves 



where 
Also define 



dtW = —iCqW + iF 



2v\w\ + vw + \w\ w . 



dcf —lur I \ 1 _ dcf 

wi = e 2 "j^iiQ — and w = w — w\ 
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SO that w satisfies 

dtw = —iCqW + iF 



where now we write F as 



I I O / \ O I I O / \ 

F = 2v\wi + w\ + v{wi + w) + \wi + w\ [wi + w) . 
Since Cq is self-adjoint with respect to (■, ■), 

dt{CqWl, Wl) = 2{CqWl, dtWl) = 2{CqWl, iCqWl) = . 



By Lemma [3.41 

^ {CqWl{t),Wl{t)) = {CqWo,Wo) < \\wo\\hi . 

Hence, there exists a constant ci > such that 

(3.5) \\wi{t)\\Hi < cih foralH. 

It can be checked by direct computation using CqV^ = iv^ that 

(3.6) P O Cq = CqO P = UJ{W,V3)V3. 

(An abstract argument using the fact that Cq preserves Vq can be given to justify the 
first equahty, which is, in fact, all we use for now). Let 

w = w — Pw 

so that Pw{t) = for all t (and hence Lemma (3.41 will be applicable to / = wit)). 
Using (13.61) . we find that 

dtW = -\iCqW + (/ - P)iF . 

Using the self-adjointness of Cq with respect to (■, ■) and the above equation for dtW, 
we find that 

dt{CqW, w) = {CqW, dtw) = 2{CqW, -\iCqW + (1 - P)F) 

= 2{CqW,F). 

Let [0, T] be a time interval over which w remains 
(3-7) < cih, 

where Ci is given in (13. 5p (so that we know a priori that w is at least no worse that 
Wl, although of course we want to show that it is better). [All future instances of < 
mean "less than a constant times" , where the constant depends upon Cq (in Lemma 
13. 4p and Ci.] This gives an estimate for F: ||F||ioo jji < h"^. By integrating, for 
< t < T, we have 

I {Cgw{t), w{t)) \ < T \\w\\L^,^^Hih^ ■ 
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Taking the sup over t E [0,T] and employing Lemma [3 ■4[ 

ll^&llroo Ml ST h'^\\w\\ roc rrl 

and hence 



(3.8) m < Th^ . 
From this, we need to infer a bound on w. We compute 

and thus, on [0,T], we have the bound 

\dM^.v,)\<h\ 

Integrating in time, we find that for all t G [0, T], 

(3.9) \uj{w{t),v^)\<hH. 
Now we perform a similar computation for uj{w,Vi). 

= W3) + {F,Vi) = luj{w,V3) + {F,Vi) . 

Appealing to (13.91) . we obtain the bound 

\dtuj{w,v4)\ <hH + h\ 
which, integrated in time, yields 

(3.10) \uj{w{t),Vi)\<t^h'^ + th^ . 
The estimates (13.91) and (13.101) give 

II^wIIl- m<h^T^ + h^T. 

II ll-'^[0,T] a: ~ 

Now, provided T < /i^^/^, we obtain 

II " [O.T] a; — 4 ' 

and thus 

l|f^||r°o rrl = ||it;||r^ m + ll-PtyllL^^ m < \cih , 

II ll-'^ICT] a: II "^[0,T]^ " " [0,T] 2: — 2 ^ ' 

for h suitably small (in terms of the constants Cq and Ci). Thus we have that the 
bootstrap assumption (13. 7p indeed remains valid over [0,T], and moreover, that 

\\w{t)\\Hi<t{l+t)h^ 

holds over the whole interval [0,T]. □ 
The following corollary is useful in streamlining Theorem [1) 
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Corollary 3.5. Suppose the hypothesis of Theorem\E holds, and that X satisfies \\ — 
A| < /i^. Then, in place of (13. 2p . we have 

l\u(t) - e^*^'/2 {v-, + e-^*^^..i/;o) < C(l + tfh' . 

Proof. It suffices to show that for any / such that -Pa,?/ = 0, we have 

(In fact, this is stronger than necessary since ||u^o||hi ^ h. The dominant error term 
arises from the fact that ||e**'^^/^fA — e**'^^/^f < h'^.) Let u{t) = e~**^^'«/ and 
u{t) = e"**^^'"?/. Henceforth we will drop the q subscript. Then, 

dt{u — u) = —iC\{u — u) + i{C\ — Cj^u . 

Note that 

(£a - C~,)u = i(A^ - \')u - 2{vl - vl)u - {vl - vf)^ . 
As in the proof of Theorem [2l we compute 

dt{C\{u ~ u),u — u) = 2{Cx{u — m), dtiu — u)) . 
Substituting the above formulas and using that iCg) = 0, we obtain 
dtiCxiu -u),u-u) = 2{Cxiu - u), (i(A2 - \^)u - 2{vl - vl)u - (vl - t;?)^)) . 

The next step is to write out the operator Cx, and address the above expression 
term by term. For the d"^ term, integrate by parts once, and then apply the Cauchy- 
Schwarz inequality. For the 6o term, use that \g{0)\ < ||5'||_h-i. For all other terms, 
directly apply the Cauchy-Schwarz inequality. The resulting bound is 

\dt{Cxiu-u),u-u)\ < |A - A|||m - M||^i||{t||ic»^^^^^i . 

Following the argument used to obtain the bound (13.51) in the proof of Theorem [21 
we conclude that 

uniformly for all t. Using that UttH/fi < \\u — u\\h], + \\u\\hi, we obtain 

\dt{Cxiu -u),u- u)\ < h'^Wu - uWjji + h'^Wu - uWhiUWhi ■ 
Integrating over [0,t] and using that n(0) = u{0), we obtain 

(3.11) \{Cxiu{t)-uit)),uit)-um<th'\\u-^l^^^^H],+th'h 
By Lemma [3.41 

II (u-u)- Px{u - u) 11^, < {Cx{{u -u)-Px{u- u)) , {{u-U)-Px{u-u))). 
By Cauchy-Schwarz, we deduce the bound 

(3.12) ||n-n||^i < (^(n - n), n - n) + (||n - n||^i + UPaC" - S) H^i) ||PA(n - S) ||^,i . 
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Note that 
(3.13) 



Px{u -u) = -Pxu = -(Pa - Px)u - P^u , 



Since P\o e = e o P-^^ 



P-,u = e-^'^-^P-J. 



This is just the evolution of the (generahzed) kernel, and hence 

\\P-,u\\Hi < t\\P'J\\Hi = t\\iP-, - Px)f\\Hi < th'Wfy^^ . 

Similarly, 

WiPx - PxY^Wm ^ h'Wuy^^ < h'{\\u - uy^ + WfWmJ ^ 

which, together with the previous bound, gives (see fl3.13p ) 

II^A(n - u)\\hi <h'\\u- u\\h. + (1 + t)h'\\f\\l^ . 

Combining this with (13.121) . we obtain the bound (absorbing /i^ilM — -ull^ri on the left 
side) 

h - ^Wli < {^X{U ~U),U-U) + {1+ t)h^U - u\\hI UWhI + (1 + tfh'WfWlr . 

Combining this with (13.111) . provided th"^ ^ 1, we obtain the bound 

h - nWl^^^^Hi ^ (1 + t)h'\\u - uWmJfWm^ + (1 + t)'/^1/||l,i , 
from which it follows that 

Mt)-m\m<il + t)h^fy^, 

which implies the estimate in the corollary. □ 

4. The free case 

We will discuss the case of g = and the initial data close to a stationary soliton 
sechx. Very precise information can in principle be obtained in this case using the 
inverse scattering method [31] - see also [6], [Z],[8]. However we are not aware of any 
reference containing that information - see [TU Appendix B] for a discussion. 



4.1. Spectral theory of the linearized operator. The explicit spectral decom- 
position of the operator Hq was discovered by Kaup [19] (see also [30] for a recent 
discussion and generalizations). We now present it in a way which will make the 
spectral decomposition of Hg natural. 
Spectral theory of operators of the form 



(4.1) 



H 



1 


'-dl + 1 " 






V2' 


2 


dl-i_ 


+ 


-V2 
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was studied systematically by Buslaev-Perelman [1] and Krieger-Schlag [21]. Despite 
the non-normality of H (if V2 7^ 0) a spectral decomposition is available once the 
existence and properties of the four dimensional set of solutions to 

is established. They are characterized by their behaviour as x — 00: 

/i2 = ^2^2, fc^O, /i^O, 

see [21], Section 5] for a careful discussion. 

For the linearization of the cubic NLS the set of four solutions can be found explic- 
itly and it is given by 

(4.2) {^+(-,A;), ^+(-,-A;), ^-(■,^/i), ^_(-,-z/i)}, k^Q, /i , 
where 

*_=ai*+, /i = (fc2 + 2)1/2 _ 

Since aiH^ai = —Hq, we have that 

{^4-,k), ^-i-.-k), ^'+(-,2/i), ^'+(-,-2/i)} 

is a basis for the solution space to HqiP = — (A;^ + l)ip. 
Now let 

V9+(x, /c) = 0-3 ^'+(x, fc) , v^-(x, /c) = criy?_(x, fc) , 
Then since a^HQa^ = Hq, we have that 

{^+(-,fc), (/?+(■,-/;;), (/?_(-, i/i), (p+{-,-ifi)} 

is a basis of the solution space to HqiP = {k"^ + \)ip. Note that 

V9„ = 0"iV9+ = (Ti(J3\E'+ = — Cr3Cri\i/4. = — (J3\E'_ 

Finally, using that aiHQiJi = —Hq, we obtain that 

{ip^{-,k), ip^{-,-k), ip+{-,ifi), (^+(-,-z/i)} 

is a basis of the solution space to Hq^P = —{k"^ + l)ip for k 0. 

The eigenvalue corresponds to a generalized eigenspace span{^i, "^2, ^3, ^4}, to 
be described now. Let ?7(x) = sechx and 

ei = —d 62 = ix 63 = -i 64 = 9x 

""^Iii Appendix |A] we show how this solution can be guessed by performing a simple numerical 
experiment even if, as we were, one is ignorant of the inverse scattering developments. We are 
grateful to Galina Perelman for explaining to us the structure of solutions to linearized operators in 
the completely integrable case. 



tp{x) ~ e 



±ikx 



(4.3) 



^+(x, k) 



(tanhx — ik) 
— sech^ X 
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Then ci (translation) and 62 (Galilean) are symplectically dual and we have fl2.2ip 
Let 

vl/, = i=-'U{e, ■ r,) , 

where U is given in fl2.13p . Then 

Ho^i = Ho^3 = 0, = ^1, Ho^4 = ^3 

dcf — 

The generalized kernel of Hq is spanned by (pj = (73'^ 

H*Qipi = HqLP3 = , HqLP2 = ipi , H^Lfi = LPs . 

Since Wa^U = 02 we also see that 

P2{xY^x{x)dx = / pI{x)'^2{x) = 1, 



Finally we recall that Hq has a simple threshold resonance given by the explicit 
formula (see Chang-Gustafson-Nakanishi-Tsai [3l §3.7]) 



(4.4) 



tanh^ X 
— sech^ X 



corresponding to k = 0. Following [T] and [211 Definition 5.18] (note a slight change 
of convention between this paper and [21]) we say that H has a resonance at 1 if there 
exists u G L°° such that Hu = u. The multiplicity of a resonance is the number of 
independent solutions with these properties. As we will recall below the maximum 
multiplicity is 2. Here we include eigenvalues as resonances: "true" resonances satisfy 
u e L°°\L\ 

This definition is equivalent to the more general definition based on the mero- 
morphic continuation of the resolvent. The potential sech^x is exponentially de- 
caying and the resolvent of Hq (the same operator without the potential term), 
Ro{z) = {Hq — z)~^, has a global meromorphic continuation to a three sheeted Rie- 
mann surface, with poles at ±1. The resolvent R{z) = {H—z)~^ can then be continued 
from the physical plane, E =^ C\ ((— cxo, —1) U {0} U (1, 00)), to a neighbourhood of S 
on that three sheeted Riemann surface. Near dzl the resolvent is meromorphic in A, 
z = ±(1 + A^). The analysis outlined in [1] (see [211 Lemma 5.2] and [211 Lemma 6.5] 
for detailed presentation) can be used to show that the definition of resonances as 
poles of the resolvent coincides with the definition above given in terms of solutions. 

Let Pc denote the symplectic orthogonal projection onto the essential spectrum, 
which we define as / — where Pd is the symplectic orthogonal projection onto the 
discrete spectral subspace Eq. The 2x2 matrix kernel Pc{x,y) of Pc is given by 
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Kaup's formula [19 



Pc{x, y) 



271 



1 



+ 

Once we know fl4.2p . this formula can also be derived by contour deformation and the 
fact that 



/ {Ho-z)-Uz = ld 



where T is any contour that encloses the spectrum of Hq - see [211 Lemma 6.8]. 
As claimed in [20] and [30| it can also be checked by an explicit calculations of the 
integral. 

We now put this into a form that is more consistent with one- dimensional scattering 
theory (see for instance [27]) and connects the basis with the basis of scattering 
solutions of P §2.5.1] and ^ §6]. 

Let 

(4.5) v+{x,k) = j^^-^^'^+{x,k) . 



Then Hqv^ = (1 + k'^)v^ and 



(4.6) 

(4.7) 
with 



v+{x, k) 



v+{x, k) 



^ikx _j_ ^^(J^'^Q ikx ^ 

T+(A;)e*'=^ 



as X 



— oo 

+ 00 



ixk 



T.{k)e 



as X — oo 



for > 



for < 



as X 



R+{k) = and T+{k) 
R^{k)=0 and T„(A;) 



iky 



iky 
iky 



iky 



Now let 



def 
- def 



H;v^ = (1 + k^)v^ 



- def 

V- = aiv 



Hov- 



Then 
(4.8) 

and 



2tt 



-(1 + e)v. 

{v+{x, k)v+{y, ky + v_{x, k)v_{y, ky) dk , 



(4.9) — [ v±{x,kyv±{x,k')dx = 6{k - k'), — [ v±{x, kyv^{x, k') dx = 
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4.2. The free linearized propagator. It follows from (14. 8 p that the propagator 
e~2**^op^ on the essential spectrum is represented by the Schwartz kernel 

271" Jr 

We will now study e~^'^^^^^ PcWq for wq appearing in Theorem [2l 
Proposition 4.1. Suppose that wq G iS(M) is real valued. Then 



+00 



(4.10) f\k) 



1 



V2^(l - i\k\) 
a{x, k) 



+00 



(1 + 2ikt{x) - kH{xy)wo{x)e-'^'' dx 



-ikx 



[l + i\k\f ' ' ' ' {l + i\k\f 
where we used the notation t{x) = tanha; and s{x) = sechx. Consequently, 

1 



(4.11) e-^''^°PcWo{0,t) = ^e^V^ [ Wo{x)dx + 0{t- 



.3/2N 



Proof. Let 



V+g{x) 



f+(x, k)g{k) dk . 



2tx Jk 

be the "inverse distorted Fourier transform," which gives 

1 



V:f{k) 



fl(x, k)f{x) dx . 



the "distorted Fourier transform," associated to the operator Hq. With this notation, 
we have 



= V+M{t)Vla^ + aiV+M{-t)Vl(j^ai, M{t)f{k) = e 



fik) 



Consequently, for wq real, 

^-litHop^^^ = [l O] P^e"^^*^ 



Wo 



[1 0] {V+M{t)V;a3 + aiV+M{-t)V;a3(Ti) 

1 



Wo 



[1 0] {V+M{t) + aiV+M{-t))V; 



Wo 
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We write the function v^{x, k) given by (14 .Sp as 



v+{x, k) 



a{x, k) 
b{x, k) 



Akx 



— s{x) ikx 



b{x, k) 



and define 
(4.12) 

Tfien, 



f{k) ( VI 



Wo (k) 



27r A. 



(l + z|A;|)2 
a(x, /c) — b{x, k))wo{x) dx . 



e-2**^Pew(0,t) 



/27r 



)/(A;) dk 



(l + z|A;|)2 



fik) dk 



-^65^^/(0) + c(r 3/2) 



by the method of stationary phase. 



□ 



4.3. Nonlinear perturbation theory in the free case. Let us take a particular 
example: 



Wo{x) 



1 + h 



sech 



X 



— sech X 



l + 2h \l + 2h^ 

(the choice of scaling was made so that sech x is selected as the nonlinear ground 
state by Theorem [2]). In this case we compute /(O) = which gives 



2 '^ + ^1^^3/2 



We also have P3W0 = lj{wo,V4^)v3 = and P4W0 = cj(wo,f3)f4 ~ 0.4 ■ h'^V4. Thus, for 
t ^ 1, we have 

e-5**^op,wo(0, t) = 0.2 ■ thH - e^'^ + 0{h?) + O 

We can now apply Theorem [2] to see that the solution of 



lUt 



-Ma;x/2 - I^I^M, M(a;, 0) 



l + h 
l + 2h 



sech 



X 



1 + 2/1 



satisfies 



(4.13) e-'"\{Q, t) = l- e-^'^ e^''h + Oit'h^) + O , 1< t < h^'^^ . 

Figure [5] compares this asymptotic expression with the numerical solution. 
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h=0.1 




10 20 30 40 50 60 70 80 90 100 



numerical solution 

stationary phase approximation 

h=0.3 



1 1 \ 1 1 1 1 1 r 




Figure 5. Breathing patterns for q = and u{x,0) = (1 + h)/{l + 
2h) sech(x/ (1 + 2h) (the rescahng is rigged so that the resulting ground 
state is simply sechx): we show |M(0,t)| and the asymptotic prediction 
for h = 0.1, h = 0.3. The agreement is remarkably good for times much 
longer than given in the theoretical result. A more precise statement 
can be obtained using the inverse scattering method. 



Remark. We should stress that a more precise result valid for all values of h can in 
principle be obtained using the inverse scattering method - see [HI Appendix B] and 
references given there. It would be very interesting to compare those exact expressions 
with our rough asymptotics. The results of [T^ Appendix B] show already that (14.131) 
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can be corrected since we know that 



where 



u[x. 



t) =e*'^('*)sechx + CLoc (^-^^ 



( sin^ vr/i 
log 1 + 







c 



cosh^TrC; C' + (l + 2/i)2 



2 u1 



(^sech ttC 



c/C ^ 0.6 



dC 







'0 1 + 

Hence, in the apphcation of Theorem [2] the error terms 0{h'^) in (14.131) are optimal. 



5. Small external delta potential 

In this section we will use Theorem[2]to prove Theorem [1] stated in the introduction. 
For that we will follow the same path as in ^ and provide a spectral decomposition 
of the linearized operator with the 5o potential. The scattering coefficients, R± and 
T, appearing in (14.61) and (14. 7p are now more singular which makes the asymptotic 
analysis more complicated. 

5.1. Basis of solutions to Hqip = ±(A;^ + 1)^. Using the Kaup basis (14. 2p for the 
free problem we find a complete set of solutions ip to the equation Hgip = {k^ + 
where 











2 sech^(a; + sgn(a;)6') 



-2 -1 
1 2 



2q 



So 






-So 



with 9 = tanh"^ q, see (l233|) . 

Let s = sech(x + sgn(x)^), t = tanh(x + sgn(x)^), and /i = (fc^ + 2)^/^ > 0. With 
unknown coefficients A{k), B{k), C{k), and D{k), we look for ip{x, k) of the form 
(5.1) 



it -iky 



it - f^f 



gM(-e) \ ^0 



B 



\t-ikf 

o2 



\t + ikf 







it + 





X, 



For the unknowns A{k), B{k), C{k), and D{k), two equations are obtained by 
requiring continuity at x = and two more equations are obtained by requiring the 
appropriate jump condition in the derivatives at x = 0. This gives rise to the 4x4 
system analysed in detail in Appendix O 

By comparing ip{x, —k) and ip{x, k) asymptotically as x ^ — oo, and noting that 
both solve Hgip = {l+k'^)ip, we find that ipix, —k) = ip{x, k) and hence A{—k) = A{k), 
and similarly for B, C, and D. 
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Here is a typical consequence of the formulas from Appendix O An eigenvalue at 
1 + k'^ comes from finding a solution to 

with ImA; < 0. More generally, a solution will give a resonance or a pole of the 
resolvent. The following lemma is derived from the computations in Appendix O 

Lemma 5.1. For g < 0, < |g| <C 1 the operator Hq has one eigenvalue, fi^ , near 
±1, 

(5.2) ^f = l-q\ 

The corresponding eigenfuctions can he chosen to he real and satisfy ctiU^ = , 
where 

\21 



Consequently, 



1 1/2 



(tanh(|x| +6) - q)' 
— sech^(|x| -(- 



g<?(kl+e)^ ^ = tanh-ig. 



(5.3) \Uq{x)\ < C|g|5e-I'?^l , j m^(x)* a^u^{x) dx = I . 

For < g <^ 1 the operator Hq has no eigenfunctions near ±1 and the thresholds ±1 
are not resonances. 

Remark. The normalization of is consistent with the spectral decomposition of 

Hq - see 

We next analyse what happens when k = —i and /i = 1. Since explicit formulae in 
that case do not play a role in our analysis, the spectrum of Hq (and equivalently of 
Fq = —iCq) near zero is analyzed by more general methods in Appendix [Bl There we 
proof the following lemma: 

Lemma 5.2. For < |g| ^ 1 the generalized kernel of Fg is given hy {ivi,dxVx\x=i}, 

Fq{ivi) = , Fg{dxvx\x=i) = ivi . 
In a neighhourhood ofO, Fg has two eigenvalues 

± ^ f ±g5 + 0(g3/2) ^ > , 

" \ ±t\q\^ + Oi\qf/^) q<0. 
The two eigenfuctions, w^, are odd, and satisfy cr^w^ = w^. 



Remark 5.3. Note that Hg = 2F„, and thus the eigenvalues of Hg occur at 2A^. 



We also see that there are no embedded eigenvalues in the continuous spectrum: 
they would correspond to real poles in A and D. Hence for /c G M \ {0} the solutions 




q>0 



no threshold res 



ess. spec. 




ess. spec. 

evals at ±2q^l'^ 



2d gen. kernel 



Figure 6. The spectrum of the operator Hq for g < (repulsive 5 
potential) and g > (attractive b potential). The threshold resonances 
of the free problem become eigenvalues for the repulsive potential which 
is counterintuitive. 



■{/^(x, /c) and '?/'(x, —h\ or the solutions '?/'(x, /c) and ip{—x, k), form a basis of tempered 
solutions to 

Hgu = {k^ + l)u. 

Our operator Hg is the Hamiltonian matrix for the quadratic form given by 



L = JHr, 





-dl + l 



-dl + 1 




+ 2 sech^(x + sgn(x)^) 



]-l 2 
2 -1 



-2q 



5o 
Sol 



but all we need are the general structural properties. 
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5.2. Spectral decomposition of Hq. Let be the symplectic projection on the 
symplectic orthogonal of the discrete spectrum of Hq - which we know consists of 4 
eigenvalues for g > and 6 eigenvalues for g < 0. 

As in the case of g = 0, we want to write the Schwartz kernel of as 



2ti 



'v^{x, k)v^{y, k)* + v_{x, k)v_{y, k)*)dk . 



where 

and 
(5.4) 



Hqv± = ±{e + l)v± , H*v± = ±{e + l)v± 



1 

2^ 



v±{x, k)*v±{x, k')dx = 5{k — k') 



1 

2^ 



v±{x, k)*v^{x, k')dx = . 



Now the generalized eigenf unctions are in fact double with ±/c corresponding to the 
single generalized eigenvalue /c^ + 1. 

A comparison with standard one dimensional scattering - see [271 (1.26), (1.30)] - 
shows that the states f+(x, k) should be chosen so that they satisfjo (14. 6 p and (14.71) 
- see [H §2.2.2] and [2T1 Proposition 6.9] for a full justification of this in the case of 
the system ( 14. ip . 

We compare these asymptotic formula to the properties of ip{x, k): 



ip{x, k) ~ 

and 

ip{—x, —k) ~ 
This shows that 
(5.5) 
where 



B{k){l - ikfe'^'' + C{k){l + ikfe-'^'' 
(1 +ifc)2e^^■^ 



X 
X 



-oo 
-oo 



(1 -^A;)2e^^^ 

B{-k)ll + 2fc) + C{-k){l - ikfe-^^"" 



X 

X 



+ 00 

— oo 



v+{x, k) 



def r a+{k)i}j{—x, —k) k > 



a^{k)i/j{x, k) 
1 



k < 



[l±ikyB{Tk) 
Note that v+{—x, —k) = v+{x, k). Define 



f_(x, k) =^ crit>+(x, k) 



and 



V±{x,k) =^ Cr3W±(x, k) . 



The states with ±fc > correspond to e± in the notation of [27] , 
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5.3. Propagator exp{—itHg/2). The continuous spectrum part of the propagator 
appearing in Theorems [T] and [2] can now be written as 

e-^''"^P,{x,y) = — [ (e-5**(^+^')w+(x,A;)i)+(i/,/t)* + e^'*(^+'='V(x,/t){}_(l/,A;)*)rfA;, 



27r 

where v± are given in § 15.21 We have the analogue of the first part of Proposition 14.11 
Since the proof is exactly the same, it is omitted. 

Proposition 5.4. Suppose that wq G 5(M \ {0}) fl L°°(R) is real valued and even. 
Then 



1, 



V 2vr 



(5^6) 2 



f{k) = —= / {a^^{x,k) -hc^{x,k))wQ{x)dx , 
V In 



where a and b are defined by 
and acv{x, k) = {a{x, k) + a(— x, k))/2, bcv{x, k) = k) + b{—x, k))/2 



In the above proposition, we reexpressed integrals over M as integrals over (0, +oo) 
in (15. 6p using that a{x,k) = a{—x, —k) and b{x,k) = b{—x, —k). This implies that 
f{k) is even, and that ac^i^x, k) and b^^i^x, k) are even in both x and k. 

5.4. Asymptotic analysis of the breathing patterns. We will now prove (11.51) 
by describing the asymptotics of 

(5.8) w{x, t) = e-5^*^«p;tt;o(x, t) , 

at X = 0. Here wq is assumed to satisfy (11.31) . In particular, 

(5.9) w{Q, t) = y?^"(a(0, A:)e-^**(^+'=') + 6(0, k)e"^''^^+''''>)f{k) dk , 

where a and b are defined by (I5.7P and / is given in (15.61) . 

We now focus on the form of the expression for f{k) and derive some of its smooth- 
ness and decay properties. The behaviour of / for large values of k can be de- 
duced directly from the definition of f{k) as the pairing of wq with a solution ip to 
Hgij = (1 + k^)^. 

Lemma 5.5. For real and even wq satisfying (11.31) . and for f defined by (15. 6p we 
have /Ik^ G C°°(]R±). For \k\ > e we have 



1 + P 
uniformly in q. 
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Proof. We recall from (14.121) (the formal structure of f{k) is the same as in the free 
case) that 



Wo {k) 



v+{x,ky 



Wo[X) 

-wo{x) 



dx . 



[l + k')f{k) 



dx 



The formulae for a(x, k) and k) above, and the formulae in Appendix [Cl show that 
the k) are uniformly bounded in x and in A;, for \k\ > e. Since {k"^ + l)v+{x, k) = 

HqV^{x, k), integration by parts (see the formula for Hg in (12.131) ) shows that 

2? .N* K(0)" 

"-dl + l-4v^ -2v^ 1 lwo{x) 



f+(x, k) 



Wo[X) 

Wo{x) 



v+iO, ky 



f+(x, ky 



dx 



+ 



v+i0,kyi2q 



,^o(0). 



+ 



K(o- 



<(o+)- 

^o(0+). 



where the last term came from the fact that wo{x) = u{x, 0) — vi{x), u{x, 0) ^ H^, so 
that wq is continuous at x = 0, and the Wq{0±) terms come from integation by parts. 

The right hand side is uniformly bounded for > e which proves the lemma for 
p = 0. 

We can now proceed by induction noting that 

Hgdlv+{x, k) = dlHqV+{x, k) 

= (/c^ + l)dlv+{x, k) + 2pkdl~^v+{x, k) + p{p - l)dl~'^v+{x, k) , 

and that for \k\ > e, \d^v+{x, k)wo{x)\ < C,, x G M. □ 

We now derive a workable expression for f{k). The formula ( 15. ip and (15. Sp show 
that for k > 0, we have 



a{x, k) 



{t - ifc)2e^'=(^'+^) A{-k)s'^e 



B{-k){i + tky B{-k){i + iky 

it - zA;)2e*^(^-^) ^ C{-k){t + ikfe^'^^''-^^ D{-k)s'^e^'^''~^^ 



and 



6(x, k) 



{l + ikf 



B{-k){i + iky 



B{-k){i + iky 

A(-fc)(t + /i)'e-^("+^) 

B{-k){i + iky 



B{-k){i + iky 



X 



xl + 



C(-A:)sV*^<"-^) ^ D(-A:)(t-/i)V("-^) 



;i + iA:)2 5(-A;)(l + 2A;)^ 



B{~kyi + iky 



x_ 
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where s = sech(a; + sgn(a;)^), t = tanh(x + sgn(a;)^), 6 = tanh~"^(g), and n = + 
From these expressions, we deduce that for x > 0, k > 0, 

(1 + C{-k)){t - ikfe'^^''+^^ {t + zA;)2e-^'=(^+^) 



(^cvi^: k) 



2B{-k){l + ikY 2{l + iky 

{A{-k) + D{-k))s^e-''^''+^ 
2B{~k){l+tky 



hAx, k) - 2B{-k){l + ikY 2(l + zA;)2^ 2B{-k){l + ikf 



and thus (using that A{—k) = A{k), etc.) for /c > 0, we have 



where 



h{k) = -j= I {{t + ikf + s2)e-*^("+')wo(x) dx , 

\/2'K 



(5.11) oo 



2ti 

By differentiation under the integral sign, integration by parts, and Taylor's theo- 
rem, we have 

Lemma 5.6. Let fj be defined by (15.111) and suppose wq satisfies (11.31) . We have, for 
each i = 0,1,2 .. ., 

Il/f lU^ <?(i + l^l), 

with the implicit constants depending only upon i and Wq ( specifically, the "q " on the 
right side could be replaced with a finite sum of seminorms of Wq). Moreover, 

fm = m + kf[{0) + eg{k) 

where g{k) is a smooth function satisfying for each £ = 0, 1, 2 . . 

ll^^'^l|L^<g(l + l^l)- 

Now we return to the computation of w{0,t) given by ( 15. 9p . Because of continuity 
at a; = we conclude that 

_ 1 / (g-^fc)V^-^ A{-k){l-q^)e~ ^^~ 
(5.12) a(U, k) - j^^, I B{-k) + B{-k) 

and 6(0, k) = bi{k) + 62(^)5 where 

e-^^ gVfe^ + A(-fc)(g + /i)V^^ 

^'■'^^ ^'^^^ = -JTTWBhk)' '^^^ {l + ^kyBi-k) 
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Upon substituting fl5.12p . (15.131) and (15.101) into (15.91) . we obtain an expression with 
many terms. We first observe in the following lemma that, fortunately, many of these 
terms are of lower order. 

Lemma 5.7. Forwo satisfying (11.31) . and fi{k), f2{k) defined in (15. lip , we have that 
each of the following 

oo POO 

e-'^'"''a{0,k)f{k)dk, / e'^'"''b2{k)f{k)dk, 

Jo 

°° i...2 &(0,fc)(A(A;) + D(fc)) 
' 2B{k){l-ikY 
is of size 



t^/v ^ \tyv ■ 

We will prove this lemma later. In the next lemma, we deduce the asymptotic form 
of the dominant terms in the expression for w{0, t). 

Lemma 5.8. For wq satisfying (11.30 . and fi{k) defined in (15. lip , we have 

e^'^^' (l + C{k) 



(5.14) 



-fi{k) + f,{-k)\dk 



/o 2(1 + P)25(-A;) V B{k) 

Combining Lemmas 15.71 15.81 we obtain the following proposition. 
Proposition 5.9. For w{0,t) given by (15. 9p we have for t ^ 1, 

Remark. The leading expression in Proposition 15.91 is formally the same as the 
expression in the case g = in §4.3[ For the case described in Fig. [H 

wo{x) = sech ( — - — ) — sech(|x| + tanh""*^ q) , 

1 + q \l + qj 

and we have 

/•oo 

wo{x) = 2q. 



Now we develop some preliminaries in order to prove Lemmas 15.71 and 15.81 To 
streamline the presentation, we introduce a definition: 

Definition 3. A function h{k,q) is conormal (at k = 0) uniformly in q, 

he A, 
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if for <i<4:, 
(5.15) 



\dih{k,q)\<C, for \k\>l, 
\ikdkYh{k,q)\<Ce for \k\ < 1 , 
with the constants independent of q. 



We note that the sum and product of conormal functions is conormaL The two 
main types of lower order terms that we encounter arise from either q times a conormal 
function or k'^ times a conormal function. The former will give an error of size q^ /t^^"^ 
and the latter an error of size q/t^/"^. This will follow (as we will see in more detail 
in the proof of Lemmas 15.71 15.81 below) from Lemma 15.61 and the following lemma 
applied with f = fj, j = 1,2 defined in (15.111) . 

Lemma 5.10. Suppose that h{k,q) is conormal in the sense of Definitionl3 Then 



(5.16) 



h{k,q)f{k) dk 



< 



3=0 



(5.17) 



6 2' 



'k^h{k,q)f{k) dk 



i=o 



with the implicit constants independent of q. 



Proof. We begin with (I5.16p . Let s = k\/t. Then the integral to be estimated takes 
the form 

^.r-'^M^)/(^)- 

Let xi.^) satisfy 

(5.18) X ^ C'^((— 1, 1)) , X is equal to 1 in a neighbourhood of s = 0. 

Clearly, 

and therefore we just need to estimate 

Using that [—is^^dsYe'^^'^'^ = e^*'^^ and two applications of integration by parts gives 

/"OO 



< iil°° ^ iij ||L° 



ds 
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Distributing the derivatives and estimating (using the factor to carry out the 
integration), we obtain the bound 

E II^^^''llL°°(|fc|>f-V2) \ / ll/^^^llL°°(|fc|>f-i/2 ) 
\0<e<2 / \0<£<2 

Now we just apply (15.151) to obtain the bound (I5.16p . 
Now we estabhsh (I5.17p . Let s = k\ft to obtain 



The remainder of the proof is similar to that above, except that we need to use 
(— 'is~^9s)^e2*'^ = 62*^ and four applications of integration by parts. □ 

We shall need the following properties of the scattering coefficients A, i?, C, and 
obtained from the more precise asymptotics in Appendix O 

Lemma 5.11 (Properties of A, 5, C, D\ 1/B{k) and C{k) / B{k) are conormal, and 
in fact 

1 k 

+ qaiik, q) + k'^a2{k, q) , 



B{k) k — iq 

+ qa'i{k, q) + k!^ai{k, q) 



B{k) k — iq 

where aj G A are conormal in the sense of Definition\^ Also, 

^ = qP,(k,q), ^ = qp,(k,q), P,eA. 
With these preliminaries out of the way, we can now prove Lemmas 15.71 and 15.81 



Proof of Lemma 5.7 We shall give the proof for 

POO 

(5.20) / e-'^'^^'aiO^k) f{k)dk. 



The other integrals in the statement of the lemma are treated similarly. By Lemma 
I5.11[ we see that in the expression (15.10p . all coefficients of /i, /2 are conormal. Also 
by Lemma [5.111 and (I5.12p . we see that 

a(0, k) = qai{k, q) + k'^a2{k, q) , aj G A. 

By the alegbra property of the conormal class, (I5.16p . (l5.17p . and Lemma 15.61 we 
obtain that (15:201) is of size 0{q^/t^/^) + 0{q/t^/^). □ 

Proof of Lemma \5.8i We write ~ to mean that the two quantities are equal with an 
error of the form q times conormal or k"^ times conormal. By Lemma 15.111 
1 k 1 + C(k) k + iq 1 
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We also take the expansion in Lemma [5.61 

fiik) = MO) + kf'iO) + eg{k), 

fi{-k) = fm-knO) + egi-k). 

Substituting the above into fl5.14p and appeahng to fl5.16p . fl5.17p and Lemma [5^6] for 
the error terms, we see that fl5.14p is equal to 



oo 



This simplifies to 

foo 7,2 /-oo 



(5.2i: 



Note that 



/•oo 7.2 poo 1,2 



oo 2 /"OO /"OO -| 



+ Jo Jo s' + l 

where, in the second term, we made the substitution s t— > Ss. Thus, 

/•oo 2 / 

In f l5.2ip . make the substitution s = t^^'^k and appeal to the above formula to obtain 

□ 

5.5. Proof of Theorem [H We will now combine Theorem [2] with the results of this 
section to proof Theorem [H We start with the following lemma 

Lemma 5.12. Suppose that wq satisfies the assumptions of TheoremUl and that 

dcf 



Ao = 1 + / wo{x)vi{x)dx 
Jr 

is the nonlinear eigenvalue specified in TheoremUl Then for the projection P\^, defined 
by fl3.1l) (with q supressed in the subscript), 

(5.22) Px,{vxo-vi-wo) = 0{q^), 
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and consequently the solution, A, to Px{vx — Vi — Wq) = satisfies 
(5.23) A - Ao = 0{q^) . 
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Proof. The definition ( I3.ip means that we need to show that uj{v\g — vi — wq, ivx^) = 
0{q^) since the other term vanishes by the reality of wq. Now, using the definition of 
Ao and the fact that Aq = 1 + we see that 



'^(^Ao -vi- Wo,ivxo) = 2(Ao - q)- j vxo{x)vi{x)dx - J Wo{x)vxo{x)dx 

= 2(Ao-g)- / vxo{x)vi{x)dx - / wo{x)vi{x)dx + Oiq"^) 



wo{x)vi{x) - {Xo - 1) J dx{vx)\x=i{x)vi{x)dx + 0{q'^) . 
The estimate (15.221) follows from 

dx{vx)\x=i{x)vi{x)dx = ^dx\\vx\\l2\x=i = 1- 



The comparison (15.231) between the exact solution and the approximate one is ob- 
tained from the implicit function theorem as in the proof of Proposition 13. 1[ □ 

The lemma shows that the assumptions of Theorem [2] are satisfied for h = Cq, 
9 = {wq is real) and 

A = Ao + 0{q^) . 
We can then apply Corollary 13.51 to obtain 

(5.24) \\u{t) - e^*^o/2 (^^^^ + e-*^^o..y;o) hi < C{1 + t)Y , 

for all < t < 

now write 



e-'^^^o.iwo = [l O] e"^ 



Wo 
Wo 



'- Wq 



Wq 
Wo 



The first conclusion of Theorem [T] given in (11.41) is immediate from (I5.24p and Propo- 
sition 15.41 once we show that 



(5.25) 



[l O] e-^'^^^o,^Pd 



Wo 

Wo 



Om{q 



3/2N 



< t < g 



-1/2 



For g > we have six contributions to the discrete spectrum, while for g < there 
are four. By Lemma [5l2] the non-zero eigenvalues in the neighbourhood of zero do not 
contribute as they are odd while Wo is even. The contribution of the zero eigenvalues is 
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0{qH) by the same arguments as in ^|3l For g < the coefficients of the eigenfuctions 
(which are uniformly bounded in H^) are estimated using Lemma [5. II by 



Cq 



1/2 



3/2 



Hence f l5.25p holds and in view of fl5.24p we have established (11 ■41) . 

To obtain (11.51) we use the above estimate and Proposition [531 The combined error 
term for ^ t ^ is 



O 



|3/2 



and that is bounded by C|g|/t^/^ for t < C"|g| 



-2/7 



Appendix A. A derivation of Kaup's basis using MATLAB 

The Kaup spectral decomposition of the linearized operator was based on the con- 
nection with the Zakharov-Shabat system and the complete integrability of the cubic 
NLS, see [19] and [30]. We rediscovered the structure of his basis of solutions through 
a numerical experiment and it might be of interest to indicate how that was done. 
The original motivation was to show that the threshold resonances for the lineariza- 
tion of the cubic nonlinear Schrodinger equation (NLS) on the line are simple which 
can be done by an explicit construction of a solution to a system of ODEs. 

The explicit resonant state of the linearized operator Hq at 1 is given by 



(A.l) 



Ml 



- sech^ x"' 
sech^ X 



To show that it is simple, we need to show that any other bounded solution is a 
multiple of Ml. 

As in standard scattering theory, the four independent solutions of Hu = u/2 can 
be characterized by their behaviour as x — > oo - see the proof of [211 Lemma 5.19]. 
In particular the resonant states can only be given as linear combinations of the two 
solutions, Ml and U2, satisfying 



(A.2) 



Ml 



+ 0(e-'"), e 



U2 



0{e 



-2x\ 



X 



-oo . 



We see that mi is given by flA.ip and Ui G Once we show that U2 ^ we will 
see that the multiplicity of the resonance is one. 

As we have already seen in §4.11 the solution U2 can be written explicitly. An 
elementary calculation con&ms that 



(A.3) 



M2 



exp( 



X) 



— sech^ X 
(tanhx + 
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This shows that 

U2 = exp 



-V2a 



i-l + V2f 



0{e 



X 



-OO 



(1 + V2y 

and, in particular, that U2 ^ L°°. 

The exact expression (1A.3P was arrived at through an attempt to produce a com- 
puter assisted proof of the fact that U2 ^ L°^. 

The solution U2{x) is obtained by solving the following Volterra integral equation 
(see [2TI (5.4)], where one should let A — + and renormalize following t— > d'^/2): 

2V2{y-x) V2{y-x) 
sinh(v^(?/ - x)) 2 smh{y/2{y - x))_ 

Solving this equation by iteration could in principle show that the solution is un- 
bounded. 

Elimination of the exponential growth in the integral equation is helpful for the 
theoretical estimates of [1] and [21] and seems essential for a succesful numerical 
scheme. Thus we consider 



U2{X 



i-c 

J X 



sedi^ yu2{y)dy . 



v{x) = exp{V2x)u2{x) 



which solves 
(A.4) 



v{x) 



+ ICvix) , ICvix 



def 



2^/2{y-x) V2{y-x) 
sinh(v^(?/ — x)) 2 smh.{\/2{y — x)) 



K{x,y) 

It is not hard to see the convergence of 

00 

(A.5) ^(a;) = 5^/C" 



K{x,y)v{y)dy , 

sech^ y exp(v^(x — y)) ■ 



n=0 



in, say C^(]R), for any k. Hence showing that v{x) 0, x ^ —00 is in principle 
possible by a numerical computation. 

We easily implement the operator }C in MATLAB. The input is an array which is a 
discretised M^-valued function on [—10, 30] with N grid points. Because of the Volterra 
structure of the equation the left limit, —10, is not important. The right cutoff, 30, 
is chosen large enough to make the effect of the potential negligible. The integrals 
are computed using the built-in trapezium rule and the errors can be estimates. 
We used = 10^ which would have to be even larger for rigorous estimates, while 
experimentally it was clearly an "overkill". 

function KT = KT(u) 

mu =sqrt (2) ; 
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[M,N]=size(u) ; 

X = linspace(-10,30,N) ; 

for j=l:N-l 

y = linspace(x(j) ,30,N-j+l) ; v = sech(y) . *sech(y) ; 
ul=u(: , [j :N]) ; 

uu(l, :)=v.*(y-x(j)) .*(2*ul(l, :)+ul(2, :)) ; 

uu(2, :)=v.*(sinh(mu*(y-x(j)))/mu) .*(ul(l, :)+2*ul(2, :)) ; 

uu(l, :)=-4*exp(mu*(x(j)-y)) .*uu(l, :) ; 

uu (2 , : ) =-4*exp (mu* (x ( j ) -y) ) . *uu (2 , : ) ; 

KT(l,j)=trapz(y,uu(l, :)) ; 

KT(2,j)=trapz(y,uu(2, :)) ; 

clear uu 

end 

KT(: ,N) = [0;0] ; 

When the numerical solution obtained using (1A.5I) with n = 10 was plotted (see 
Figl?]) we noticed that the plot of the first component looked remarkably like a plot of 
—a sech^ x, a > and the fit based on the minimum of first component (experimental 
—a) was almost exact. From the operator Hq it is clear that having one component 
of the solution we obtain the other and that quickly led to the exact solution (lA.Sp . 
This then suggests the form of general solution for other values of k and /i as given 
in mi 

It would be difficult in general, and by our method in particular, to show the 
existence of a resonance. 



Appendix B. Perturbation of eigenvalues at zero energy 

Here we present the perturbation theory for Hg — z a.t z = 0. Even though we 
could in principle obtain the same results from careful analysis of the matrix A{k, q) 
described in detail in Appendix [Cl the method used here is more general and does 
not depend on explicit formulae. It is of course close to the similar study in the semi- 
classical case, see [11] and references given there. However, since the delta function is 
clearly different from a slowly varying potential with a nondegenerate minimum we 
give a selfcontained argument. 



B.l. Grushin problem. We recall that the linearized operator acting on 

Rew ^ 
Imw 



is given by 
(B.l) 



c 



def 



' 1" 




'L,+ " 


-1 




. V. 



1(1 
1(1 



dl - Qv^ - 2q6o) 
dl - 2v^ - 2q6o) 
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• 1 St component of the numerical solution 

♦ 2nd component of the numerical solution 
The exact solution 
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Figure 7. The plots of components of exp(^/2x)u2{x), following 
the numerical computation and the exact solutions. We see that 
lim3;^_oo exp {^/2x)ul{x) = (3 - 2v/2)/(3 + 2^2) ^ 0, where U2 is the 
second component of the solution. We used = 10^ grid point and 
the plot shows the sampling of 100 points. 

where v is the nonlinear ground state. We take elements of H^{M.; C) and write them 
as column vectors of real and imaginary parts giving an identification 

H^{R; C) ~ H^{R; M) © H^{R; M) 

The elements Cji] take the 2-vector form 



-V 




'0" 




"0" 




rj + XT]' 


_ 


, e2?7 = 




, 6377 = 




, 6477 = 





XT] 


V. 



The symplectic form, in vector notation, becomes 



(B.2) 



( 


Ml 




^'1 








5 







-UIV2 + VIU2) 



44 J. HOLMER AND M. ZWORSKI 

In the matrix notation, the relations (12.211) become 



LJ 




—Tj' 




LJ 




'0 " 




—Tj' 










= 0, 


-L+ 




XT] 




_ 



LJ 




0" 


= 0, 


LJ 




XT]' + 7] 




0" 







Jl. 












Jl. 



To perform spectral analysis, we complexify the space and work on 



H^(R; C) 



The symplectic form u (lB.2p extends to H, by analytic continuation (with exactly 
the same expression as in ( 1B.2I) : we do not insert any complex conjugations). 

Following the standard procedure (see [25]) we build an invertible matrix in block 
form 



Fg-z 



R 



+ 







with suitably chosen 



i?_ : 



H, R.:H-*C' 



We will select R-, R+ to be constant (independent of q and z) operators such that Q 
is invertible with inverse represented in block form as 



E 



E_ 

The components depend on q and z and have the following mapping properties: 

E_ : iJ ^ C^ 
E^^ : . 



E:H^H, 



To find i?_ and i?+ and to compute E^[z), E^{z), E^{z) and E^j^{z), we first consider 
(Fo — z)\Q.ri, that is Fo — 2; acting on the generalized kernel. Ordering the basis using 
Cj ■ r], j = 1, 2, 3, 4, we see that 



-z 1 
-z 
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—z 


1 
























1 








—z 






1 







1 









z 









—z 


1 























1 









—z 





1 








1 








z 


1 

















z 


1 


















1 

















z 


1 




z^ 



gives us R± for which Qq is invertible: 



Cl 62?? + (264^7, R+U 



P3U 



Pu = ^PjuejT]. 



This tells us that 
El 

or, more explicitly, 
(B.3) 

We also find that 



zCi 

C2 



CieiT] + zCie2r] + (263?? + <2e4?7 







C2. 





7]' z{ri + XT]') 
ZXT] 7] 









C2. 



E 



ai 
"2 

0(4 

^ is expressed as 



zai + a2 

Za^ + Q!4 







V 





—z J UXT] — J vrj' 



or in other words E^ : H ^ C 

(B.4) E° 

We use the following formula to compute E'L_^_{z): 

E% = E\ - E'^{F, - Fo)El + 0{q') 

By the Schur complement formula. Eg — z is invertible if and only if E'L_^_{z) is invert- 
ible, so we want to find z (in terms of q) such that det(-Ei_,_(z)) = 0. We know that 
: C2 ^ C2 is 



E" 



^2 0' 
z^ 



and thus have all the ingredients to analyze the perturbation. 
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B.2. Substitutions. We have 



sech^(|x| + q) = sech^x — 2gsech^xtanh |x| + ■ ■ • 

and thus (to first order in q) 

— L° = Ggsech^ xtanh |x| — q5o{x) 

Li — L'^^ = 2gsech^ xtanh |x| — qSo{x) 

and therefore (to first order in q), 

p ^ 2gsech^ xtanh |x| — g5o(2^) 

^ ° — 6g sech^ a; tanh |x| + g5o(a^) 

We will use the notation rj = sechx and a = tanhx. Using ( 1B.3I) we see that 
{Fg - Fo)El :£? takes the form 



2qzxrfa sgn x 2qrfa sgn x — qSo 

—Qqrfa^ sgn x —Qqzrfa{l — xa) sgn x + qz5Q 



From this, and (]B.4[) . we compute E^_{Fq - Fq)E^ : ^ takes the form 
E'^{F,-Fo)El = 



-qz'^a — qP 

qi + qz^^ 



where 

and thus 
(B.5) 



P = 6 7] a sgn X = 1 , ^ = 1 — 2 r] a sgn x = 



0{q') ■ 



'l + qa)z'^ + q 

(1 - q5)z'^ 



By expanding deti?^_,_(x) as see that Elj^ fails to be invertible when z = iiq^^"^ + 
0(g^/^). The explicit generalized kernel of Fg given in the beginning of [|3] shows that 
the double eigenvalue at persists under perturbation. We can now give 

Proof of Lemma \5.^ We only need to check the properties of w^. The equation 
(Tsiu^ = follows from the fact that a^F^a^ = —Fg. Since the eigenfuctions are 
simple and Fg commutes with u{x) u{—x), and because of the as symmetry, they 
are either both odd or both even. Schur's formula (see for instance [251 §1]) shows 
that 

(B.6) Res,=,±(F, - z)-' = Res,=,± Eliz)El4z)-'El{z) , 

and we note that 



(B.7) 



Image Res^^;^±(Fg - z) ^ 
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Eliz) 



'Ci 




—f]' Z^t] + XT]') 




'Ci 


C2. 




ZXT] 7] 




C2. 



Ml 
U2 



—z f Uixr] — f 

Z J U2{xr])' - J UiT] 



+ OH{q\Cy) 

(g 



This, ( iRHjl . and (iRTl) show that 



7] 

iiq^^'^XT] 



Oh (q) 



Hence is approximately odd, and consequently odd. 



□ 



Appendix C. The system of equations for A, B, C, D 

Here we describe how to solve for the coefficents A{k), B{k), C{k), and D{k) in 
(CT . Define 



fix,k) 
g{x,k)^ 



dcf 



■ip{x, k) 



Set A = e^^^-^^^M, B = e^'^^B, (7 = C, and D = e^^^-^^^D, 9 = tanh"^ q. Denote 
/(0±) = lima-^oit /(a;), etc. Using that s(0) = (1 — g^)^/^ and t(±0) = ±g, we obtain 

e^'V(0-) = (g + ^fc)'-i(l-g') 

e*^-7(0+) = B{q - ikf + C{q + zA:)^ - D(l - q^) 

e^^Vo-) = -(i-g') + i(g + /i)' 

e^'=^^(0+) = -5(1 - g2) _ ^(1 - g2) + + 
The two equations we obtain by requiring continuity at x = are 



(C.l) /(O) "=^^(0-) = /(0+), and gi^)"^ g{Q-) = g{Q+) 



jke fi 



/'(o- 
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We further compute, from the formula for ip, that 

= (g + ikfik - 2(1 - g2)(g ^ ^k) + a{ - {I - q'^)fi - 2(1 - q^)q 

= B{{q- ikfik + 2(1 - q^){q - ik) 
+ c(^-{q + ikfik + 2(1 - g2)(g + ik) 
+ ^(Ml-g') + 2(l-g2)g 
= ( - 2(1 - q^)q - (1 - q^)ik) + A^- 2{q + - q') + fi{q + fi)' 
= 5(2(1 - q^)q - (1 - g2),fcj + ^(^2(1 - q^)q + (1 - g2)^fc^ 
+ D(2{l-q'){q + fi)~fiiq + fif^ 
The form of the derivative compatibihty conditions is: 
(C.2) /'(0-) - /'(0+) = 2g/(0) , g'{0-) - g'{0+) = 2qg{0) . 

The four equations ( 1C.1[) .( |CI2!) give rise to the 4x4 system 



Jke I 



Jke t 



9(0- 



(C.3) 



A{k,q) 



' A ' 




Uk 


B - 1 


= Q 





C 


-2(^ - q) 


D 




-2(1 -g^). 



with the coefficient matrix A{k, q) given by 



l-g2 

(1 - q') 



(g — iky 



{q + ikf 
l-g2 



{q-ik){^-q) {q + ik){fi-q) 



(1 - q') 



-(g + /i)(A;2 + g2) (I - q^)(^q - ik) (l-q^)(q + ik) -(q + ^)(k^ + q^) 

C.l. Exact solutions. The solution is obtained from Mathematica or, in principle, 
by Gaussian elimination. Recalling that 9 = tanh~^ g, we have 



A 
B 
C 
D 



2e^-^^+i')\{ik + g) (-1 + g2 



1 + g2 (-2 + A;2 + 2g2) + 2g (P + g^) /i + (F + g2) ^2 
^-2ikeQ^ _ -q^i^- ^ ^(-g + /i) - ig(2g + /i))(fc(g + /x) - ^(1 + g/x)) 
(1 + g2 (-2 + P + 2g2) + 2g (P + g2) /i + (P + g2) ^2) 
^g (-1 + (2 + P) g2 + 2g (P + g2) + (p + ^2) ^2) 
A; (1 + g2 (-2 + A;2 + 2g2) + 2g (P + g2) /i + (P + g2) 

2e(-*^'+^)^g(zA; + g) (-1 + g2) 
1 + g2 (-2 + P + 2g2) + 2g (P + g2) /x + (P + g2) 
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The numerator in the expression for B is {k — iq)v{k)w{k), where 

v{k) = {i + kjj) + q{k — ijj) — 2iq^, w{k) = {—i + kjj) + q{k — ifj) . 
We clearly see that k = iq is a. root of B, and we further find that aX k = iq, 
A{iq) = 0, B{iq) = 0, C{iq) = 1, D{iq) = . 



Thus, 



solves the equation 



(tanh(|x| +6) - qf 
-sech^(|x| + 6*) 



giving an eigenvalue when g < 0. 

We will now specify a branch of /i = v2 + P^, and study the roots of v{k) and w{k) 
to check for consistency with Appendix [Bl Since 2 + k'^ has roots at ±^^/2, we will 
cut along the imaginary axis, and take fi as the branch defined on the domain 



(L\ —ZOO, — +ioo) 

that is real and positive for A; > 0. 

We now examine v{k) for < |g| <C 1. Setting k = 
11=1- inq^/"^ + K^q + Substituting yields 

v{k) = -2i{k^ + l)q + 0{q^'^) 



-i + Kg^/^, we find that 



and thus a root occurs at k, = ±i, i.e. when k = —i ± iq^^"^ + 0{q). Substituting 
k = —iiiiq^^'^ into the numerator of the formula for B, we obtain (9(g'^/^), while substi- 
tuting into the denominator, we obtain 0{q), and thus we have found an approximate 
root of B. This implies that we have eigenvalues at 1 + k'^ = ±2g^/^ + (9(g). The roots 
of w{k) occur near k = +i, giving nonphysical poles of the resolvent {Hg — [k"^ + 1))^^. 
From the above formulas, we have 



A _ 2ze(*'=+")^A;g(-l + g2) 

B~ {i + k{q + fi)- iq{2q + /i))(fc(g + /i) - i{l + g/^)) 

l_ _ e^^fcgfc (1 + q^ (-2 + k^ + 2q^) + 2q + g^) /x + (fc^ + q^) fi^) 

B ~ {k- iq){i + k{q + fi) - iq{2q + /i))(/c(g + yu) - i{l + g/i)) 

C _ le'^'^'q (-1 + (2 + P) g^ + 2g (P + g") ^l + {k^ + q^) /i^) 

B {k - iq){i + k{q + fi) - iq{2q + fi)){k{q + fx) - i{l + q^i)) 

D _ 2ie(^'=+")^A;g (-1 + g^) 

5 ~ {i + k{q + /i) - zg(2g + ^J)){k{q + /i) - z(l + g/i)) 
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C.2. Behaviour for large k. The behaviour for large values of k could be deduced 
from general principles of scattering theory. Here we proceed directly using the matrix 
A{k, q) which we write as A{k, q) = ^o(^) + Q^ik, q), where 

1 -k^ -1 



Ao{k) 



12^ 1 1 -12^ 

1 —ikfi ikfi 1 
— /i/c^ —ik ik —f^k^_ 



and 



■^0 





" 1 


P 


1 


-/i " 


1 




1 


ikfi 


i/k 


2(1 + 




1 


—ikfi 


-i/k 




.-1 


-k^ 


1 





fi = v2 + P^. For > e > 0, we have 

with the implicit constant in the first estimate dependent on e. Hence 

qA,'B = Oc^^c4q)- 

For q small enough, depending of e, we can used the Neumann series inversion of 
/ + qAQ^B to obtain, and consequently, for \k\ > e. 



+ OiqV{k)). 





" A ' 




' iik' 




" 4ik ' 


(C.4) 


B- 1 
C 
D 


= q{I + qA,'B)-'A,' 




-2i2 
-6 


= ?A"' 




-2i2 
-6 



(C.5) 



■ A ■ 

B-1 
C 
D 



k 



2^2 



2ik + 2/i 
-3i/k — ?>ikjj? 
"iijk — ikfi"^ 
-2ik + 2/i 



OiqVik)) 



This provides the estimates needed in Lemma [5.51 
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